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1. 


INTRODUCTION 


Locating the source of transmission is one of the tasks of military 
intelligence. In VHF frequencies and above, this is done by measuring 
the direction of arrival of the radio waves (DF) by two, three, or 
more stations, and then pinpointing the source by drawing straight lines 
on a map, 


i In HF (2-30 MHz) the situation is different since the main mode 

' of propagation is by reflection from the ionosphere. This enables the 
' use of a different technique where the azimuth and elevation of the 
incoming wave are measured together with the state of the ionosphere 
at the time of reception. 


eo All this information is presented as an input to an ionospheric ray 
' tracing program that calculates the ray's trajectory from the reception 
i site to the transmitting point. 


The actual system needed is quite complicated and elaborate. Apart 
from the equipment to measure the direction of the arrival of the 
incoming wave, there is a requirement to measure in real-time the 
structure of the ionosphere at the area of interest. It is required, { 
therefore, to measure the electron density profile and the tilts of the i 
ionospheric layers. This can be accomplished by using a Digisonde | 
together with a special "Drift Attachment" equipment that can measure 
ionospheric tilts by Doppler techniques. 


At the Propagation Research Team, Signal Processing Division, 
Communications Systems Center, CORADCOM, an attempt is made to implement 
Such a system using the existing equipment at the ionospheric station 
at Naval Ammunition Depot Earle. 


At this stage, integration of the Digisonde, Drift Attachment and 
PDP-11/40 mini-computer is carried on with the intention of adding the 
DF equipment at a later date; thus actual field measurements in real- 
time will be possible in order to determine the overall accuracy of the 
system. 


This report covers one important step in the above mentioned 
process - the computation. 


- 


A three dimensional ionospheric ray tracing program was chosen. (1) 
This FORTRAN program was originally written for a CDC-3800 machine and 
had to be adapted to the PDP-11/40 mini-computer. 


| 

| 

} 
The report is essentially an instruction book that shows how to run 

the ray tracing program on the PDP-11/40, and is divided into two main 

parts. The first part describes generally how to write into the machine, 

edit, debug, and run a FORTRAN program. The second part depicts the 

operation of the ray tracing program itself. 


How to operate the PDP-11/40 computer 4§ described in full in the 
various manuals of the machine, especially a (3 However, a reader 

might easily get discouraged because of the volume of material that has to 

be covered. The first part of this report is an easy,step-by-step 

procedure. By adhering to it, the reader will be able to operate the 

PDP-11/40 machine safely and run FORTRAN programs on it. However, it is 

strongly suggested that the machine manuals be read because the * 
information included in this report is very limited for obvious reasons. 
The second part describes the ray tracing program that exists in the 
PDP-11/40. However, it is assumed that the reader is already familiar 
with the ray tracing program. Therefore, it is necessary to read ref. 1 
in advance, especially pp 1~50, which generally describes the ; 
theoretical basis and actual structures of the program and pp 161-185, | 
which give the sample case. 


ee 


| 
| 
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2. Operation of the PpP-11/40 Machine 


The following instructions are relevant only to our machine and are 
not intended to be taken as general rules for other POP-11/40 machines. 


Our machine has only one disk (DK) and on Input-Output (1/0) device - 
the teletype (TT). ; : 


‘ 


2.1 Turning On 


Make yourself familiar with tne various switches of the machine. 
Actually, there are three panels to be concerned with: 


a. The teletype (TT panel) 
b. The main panel of the machine 
c. The disk (DK) panel 


- Locate the three switches (SW) on the left side of the TT panel. 
Make sure thac POWER SW is pressed on ON, the middle SW on LINE and 
BAUD RATE is pressed on 300. 


- Make sure that ENABLE/HALT SW on the DK panel is pressed to 
HALT and that the POWER SW (the one with the key) is in OFF position. 


- Verify that the RUN/LOAD SW on the DK panel is pressed to LOAD, and 
that the WIPROT SW is pressed the same way. 


- Turn the POWER SW on the main panel to POWER, lights will turn on, 
and you will hear the ventilators start working. After a few seconds 
you will hear a mechanical click from the DK section. This will also 
be signified (in case you miss the click) by the turning on of the LOAD 
light in the DK panel. 


- Turn the LOAD/RUN SW on the DK panel to RUN. The DK will start to 
revolve. Wait some time until RDY (ready) light on this panel turns 
on. This signifies that the DK has reached the required RPM. 


2.2 Boot-Strap 


Now the computer is ready to be activated. The action of interfacing 
the machine together with the operating system that resides in the DK 
and thus allows communication with the computer via the TT - is called 
Boot-Strapping. The action itself involves writing into the computer's 
memory a short program in machine language. The program is shown in 
Table l. : 


TABLE 1: 


Address 


901000 
001002 
001004 
001006 
001010 
001012 
001014 
001016 
001020 


Boot-Strap Program 
(All No. are octal) 


Content 


012700 
177406 
012710 
177400 
012740 
000005 
105710 
100376 
005007 
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Entering the program is performed through the various switches on 
the main panel. On the lower left side are located 18 SW numbered 0-17 
which are divided into groups of three colored by red and purple. Each 
group of three stands for one octal number which is entered in binary form. 


When a SW (one of the eighteen) is in the down position this means 
binary "0", and when in the up position-binary "1", The position of 
the groups of three SW corresponded to the position of the octal numbers 
in Table 1. So, if you want to represent the octal No. 012700, you start 
from the left (for example): 


Take the first group from the left and press the SW down; this repre~ 
sents the octal number 0. Proceeding, you take the next group of three 
SW (nos. 12-14). In this group, SW number 12 has to be in the up 
position and the other two down. 


In the next group (SW 9-11) only one SW has to be in the up position, 
SW number 10. This represents the octal number 2. The next number is 7, 
and is represented by all three SW (nos. 6-8) in the up position, etc. 
The actual Boot-Strapping is done as follows: 


- Write the octal number 0010000. 


- Activate the LOAD ADRS (Load Address) SW (on the main panel) and 
you will see this number appearing on the ADDRESS lights above the 
eighteen SW (the position of the lights corresponds to the position of 
the SW). 


- Write the content of this address (from Table 1), i.e., 012700. 


- Activate the DEP SW (on the main panel) and you will see this 
number appearing on the DATA lights just below the ADDRESS lights. 


- From now on, you will be advanced automatically. So all you have 
to do is to write the content of Table 1 and activate the DEP SW each 
time. After each step, check the content by looking at the lights. If 
you want to examine a former address, write that address and activate the 
EXAM SW. If the content is in error, load the right one. In any case, 
whether the content is right or wrong, do not forget to return to your 
original address. 


- After completion of writing Table 1, write again the address 
001000 and load it by the LOAD ADRS SW. 


~ Move ENABLE/HALT SW to ENABLE 
~ Activate the START SW. 


You will see lights flashing on the ADDRESS and DATA lights line and 
after several seconds the TT will print. 
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The dot (.) in the left corner of the page signifies that the 
computer is activated and is now ready to receive instructions, run 
programs, etc. 


In case you do not get this printing, try Boot-Strapping again 
because it is easy to err. If you are sure of your Boot-Strapping 
procedure and you still do not get the desired results, it means that 
something is probably wrong with the computer and expert help is needed. 


After Boot-Strapping is completed, it is advisable (but not necessary) 
to turn the POWER SW on the main panel to LOCK position. This 
deactivates all the switches on this panel so that accidental change 
in position of one of them will not harm the operation of computer. 


2.3 Turning Off 


The turning off procedure is important in order to avoid physical 
damage to the DK and/or loss of files on the DK. 


- If the POWER SW on the main panel is on LOCK position turn it to 
POWER after verifying that all the SW are in their proper positioning, 
i.e., that they are in the same position as after Boot~Strapping. 


- Verify that the last response of the computer on the TT is (.). 
This means that the computer is in the KEYBOARD MONITOR (KMON) mode of 
operation. If not, read further on how to transfer to this mode. 


- Change the ENABLE/HALT SW on the main panel to HALT. 

-Turn the LOAD/RUN SW of the DK panel to LOAD. The RDY light will 
go off. After about 15-20 seconds, a mechanical click will be heard and 
the LOAD light on this panel will be on. 


- Shut off the machine by turning the POWER SW in the main panel to 
OFF. 


2.4 Qperation and Communication via Teletype 


At this stage, the communication with the computer is done through the 
TT keyboard. 


2.4.1 Modes of Operation-General 


After Boot-Strapping the computer is in KEYBOARD MONITOR (KMON) mode. 
This mode is signified by the fact that after each operation the computer 
responds with a (.) at the left side of the line. In this mode, running 
programs are possible. 


Other modes are FORTRAN, EDIT, LINK, PIP (only modes concerning our 
purpose were mentioned here - see the computer manuals for further ; 
. details). ; 


In the EDIT mode, you can type in programs, edit them, and store them 
on the DK. | 


In the FORTRAN mode, you can compile FORTRAN programs. 


In the LINK mode, compiled FORTRAN programs are translated to machine 
language. 


The #IP stands for Peripheral Interchange Program. In this mode 
information can be transferred between the various peripheral devices. 


Entering into a mode of operation is done from the KMON mode by 
typing the appropriate statement. For example: To enter the LINK mode 
type 


-R LINK (CR> 


The (.) is already there to signify that the computer is in KMON 
and awaiting further instructions. 


' You type R (for RUN), leave one space, and then type LINK followed 

i by <CR>. The <CR> stands for Carriage Return and this tells the computer 
| that the statement is completed and that it is its turn to act. This is 
done by hitting the RETURN key. The carriage will return to the 


* beginning of the line. No character will be printed. 


In simple language you told the computer to run the program LINK. 

ee. This program is a part of the operating system and is already on the DK. 

| Do not forget to leave the required space. Whenever you type wrong 

' statements, the computer types an appropriate message that tells you to 
correct the last statement. 


However, if everything is right then the computer responds with a 
star (*). The whole procedure looks like: 


N 
R LINK ger > 


* 


The (*) signifies in this case that the computer is in the LINK 
mode and it awaits further instructions. In order to return to KMON 
you have to press two keys together, the CNTRL key and the letter C. 

This echoes as “\C . 
The whole procedure looks like 
-R LINK <CR 7 

* SNC 

And now you are back in KMON. 

In a similar manner, you can type R FORTRN, R EDIT or R PIP and be 
in FORTRAN, EDIT or PIP mode respectively. Returning back to KMON form 
FORTRAN, LINK or PIP modes is done by the7\C. Returning from EDIT mode 
is done differently and will be explained later. 


In order to transfer from one mode to another, you have to go through 
KMON mode. This cannot be done directly. 


2.4.2 File Names 


The general structure of a file name is fil~-nam.ext 
where: fil-nam stands for file name - a string of up to six characters 
(letters and numbers only) beginning with a letter. 


ext - stands for extension - a string of 3 characters that is either 
given by you or by the computer by default. 
, 


Suppose you want to create a file by name TAKE]. For this 
computer you have to add an extension like TAKE1.ABC or TAKE1.003. 
However, if this file contains a FORTRAN program, the extension has 
to be FOR (in order to comply with the instructions that you will be 
given later on; for different instructions see computer manuals) - like : 
TAKE]. FOR. 
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Extensions given by the computer by default that may concern you are: 


OBJ which is given after compilation of a FORTRAN program. 

SAV which is given after the link procedure (LINK mode). 

BAK which is given if you create a back-up version of your file. 
2.4.3 PIP Mode 


Again, PIP stands for Peripheral Interchange Program. In this mode 
transferring information between the various peripheral equipment of 
the PDP-11/40 system is made possible. 


In our particular case, since we have only one DK and on TT, the 
program uses are somewhat limited although important. 


In order to transfer to the PIP mode you have to be in KMON mode and 
type: 


.R PIP <cr 7 


* 


The computer responds with the (*) which signifies that it is in the 
PIP mode and awaiting further instructions. 


2.4.3.1 Viewing List of Files on DK 


In order to get a list of all the files that are on the DK, type 
*/L <CR> 


The (*) was already here signifying that the computer is ready. The (/) 
Means all, and the (L) stands for List. So what you actually said is 
List All. The <CR>is done by hitting the RETURN key, and this is a sign 
that your instruction is completed and it is the computer's turn to act. 


After this instruction is typed, a complete list of existing files on 
DK will be typed. When the typing has been completed, a (*) will be typed 
again to signify the fact that the computer is ready for further 
instructions. 


2.4.3.2 Wiewing Content of a File 


Suppose you want to examine the content of a file named REACH.FOR, 
then type: 


*TT:=REACH.FOR (CR > 


(You should type the complete name including the extension) and the 
{TT will print the content of the file. The action will be ended by (*) 
as usual. 


One important point to remember, however, is that not all the files . 
on the DK are written in the same code. Some of them are in various 
versions of machine language and trying to print them will yield 
nonsense. Wa 


In particular, all the files having the extension.SAV or .OBJ that 
were given by the computer after compiling or linking will not be 
comprehensible. 


Only the so-called Source Programs will be comprehensible. So 
whatever,the file name is, type in its name according to the way shown 
above and its content will be printed. 


If you want to stop the printing while it is still in progress type: 
on Oo 


This is done by pressing simultaneously the CNTRL key and the letter 
0 key. The printing will stop immediately and the (*) will be printed. 


2.4.3.3 Deleting a File on DK 


In order to delete a file type: 


tanne.a---/D CR) 


Which means Delete All of the file whose name is signified here by 
-~--,----, After the instruction has been carried out the (*) will 
appear. 
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For example, deletion of the file named RACK.SAV will appear on the 
Paper as: 


*RACK.SAV/D 


* 


Note that although you hit the <cr> , this action does not echo on the 
paper. The second (*) means that the action is completed. 


2.4.4 Edit Mode 


As mentioned before, creating and modifying files on the DK in this 
mode is possible. Entering this mode is made via the KNOM by typing: 


-R EDIT CCR 7 


* 


The (*) is the response of the computer. 


In this mode there is a special character $ (ESC/SEL). This is the upper 
left key and it echoes as $. This is not equivalent to the usual dollar 
sign which is the upper case (shift) of the 4 key, although they appear 
to be the same in printing. 


The role of the $ (ESC/SEL) is dual - firstly, it serves as a delimiter 
between individual instructions in instruction chains, and secondly, $$ 
in this mode means what <CR does in other modes, i.e., the instruction 
is over and it is the computer’s turn to act. After execution of the 
instruction the computer responds with (*). 


2.4.4.1 Creating a File on DK 


Suppose you want to create a file on the DK. Usually you are in KNOM, 


so in order to create the file you have to type the following sequence: 
-R EDIT<CR> 


*EWDK: fil-nam.ext$$ 
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hy ret Se type in text 


*EX$$ 


Explanation 


The first step is transferring to EDIT mode. The computer responds 
by (*) and you type EWDK which means Edit Write on DK, a file having 
the name fil-nam.ext. The $$ means execute. The computer responds with 
a (*) and you type I (for Insert) and the text. Now remember: Type the 
text as if you are typing on a regular computer card. If you have a 
FORTRAN program that you want to type in, be careful to follow the usual 
FORTRAN coding form, i.e., start with the seventh colum (leave 6 spaces), 
do not pass the 72nd column, the 6th column is used as a continuation 


mark, etc. In this machine, however u do not need line numbers (which 
are automatically given by the compiler) so do not type them in. 

After each line type in <cry . In this mode it is just another 
character (but in FORTRAN this means end of computer card). 


If in the process of typing you make a mistake and notice it 
immediately, you can correct the mistake by using the DELETE key. Hitting 
this key once deletes tye last character printed. Hitting this key twice 
deletes the last two characters printed, etc. So if the characters are 
not too far away you can delete them and retype. The whole procedure is 
shown in the following example: 


ABC=SQRT (X)*SIM(ALP PLA (MAN (ALPHA) 


If after reaching the letter P you notice a mistake has just been made, 
you hit the DELETE key 5 times and the deleted characters are printed in 
the order of appearance from right to left with the special character 
, i.e., PLA(M. When you start to retype, the second (\) appears to 
you will know that everything between the two (\)'s has been deleted. 

The whole new line is now: 


ABC=SQRT(X )*SIN (ALPHA) 


EEE 


The number of times that you can use the DELETE key is unlimited 
You can even go back to the previous line - this is possible because ] 
the whole program is just a long string of characters and the <CR? is | 
the same as any other in this mode. However, going too far back is not 
worth the effort, so just keep typing and correct the error later by = | 
using methods that will be explained further ahead. After you type the | @ 
whole program, type $$ which tells the computer to store the | 4 
information on the DK. When this is finished the computer prints (*) and | 
you type EX$$. This stands for Exit. This is the way to leave the EDIT 
mode and return to KMON. 


If however, you want to check yourself then do not type the EX$$ so 
you are still in the EDIT mode, and you can check and correct your text 
(using the ways shown later on). After corrections have been made, you 
then exit. 


2.4.4.2 Calling Sequence 


Whenever you want to modify a program that exists as a file on the DK, 
you have to be in the EDIT mode, call the file from the DK to the memory, 
modify it and store it again on the DK. Suppose you have a file named 
DEMO.FOR that you want to modify, then type: 


-R EDIT (CR 
*ERDEMO.FORSEWDEMO . FORSRSB$$ 


¢ make changes 


*EX$$ 


First you enter the EDIT mode. Then you instruct the computer to 
Edit and Read a file named DEMO.FOR from the DK (by default since we have 
only one DK) and Edit and Write on the same DK a file by the same name. 
Then you make the necessary changes according to the ways that will be 
explained later on, and you Exit to KMON. Using EX$$ again sssures that 
all the changes will be saved. After exit you are agairi in KMON which is 
signified by the (.). Another way is by typing: 


*EBDEMO.FORSRS$ 
i make changes 


*EX$$ 


This also creates a back up version that by default is getting the 
ext.BAK, i.e., DEMO.BAK 


Sometimes it is useful to have a back-up. 


2.4.4.3 Copying a File 
By the following sequence: 


-R EDIT CR 
*ERDK: GEO.007$$ 
*EWDK:ABE.FORSS$ 


*EXS$ 

You have entered the EDIT mode, you Edited and Read from the DK (only 
one in our system) a file named GEO.007 and copied it (Edited and Wrote) 
on the same DK under the name ABE.FOR, and returned to KMON. The file 
GEO.007 was not affected by the act of copying. This of course is just 
an example, but the general rule is clear. If you want to make corrections 
to the new file, first you have to exit via EX$$ and then enter the EDIT 
mode again. 


2.4.4.4 Controlling the Pointer's Position 


Before you modify your text, you have to be sure where the intended 
modification is to take place. For this purpose there is a pointer. 
The pointer does nothing but point to the exact place in the string 
of characters that comprises your program, at that particular moment. 
For example, whenever you call up a file for editing (by the calling 
sequence 2.4.4.3), the pointer is initially located at the beginning 
of the program. Whenever you are not sure where you are, or you want 
to return to the beginning, type: 


*BS$ 
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The (*) was already there as a result of previous operations. The B 
means Back and the double ESC/SEL means, in this case, execute. 


The following is a list of the most important instructions needed. 
For further details see the computer manuals. 


*+nA$$ - Advance the pointer +n lines from the current position of 
the pointer. If the pointer is in the middle of the current line it will 
nevertheless be positioned at the beginning of the appropriate line as 


though the pointer was originally at the beginning of the current line. 


As mentioned before, each line is terminated by hitting the RETURN 
key. This actually prints two instructions into the memory. First: 
Carriage Return KcR> and then Line Feed < LF» and not just <cCR7as was 
mentioned earlier, for the sake of simplicity. Now, when modifying the 
text, care should be taken not to lose any of these instructions. 


For example, when you type 2A$$ you mean that you want to advance the 
pointer by two lines. The pointer is then moved until it encounters the 
second string of <cR> <LF) and then it is stopped. This is of course the 
beginning of the next line. 


You can advance forward or backward (by the (-) sign). 
*/A$$ - moves the pointer to the end of the text. 


*OAS$ - moves the pointer to the beginning of the current line. 
(equivalent to OJ$$). 


*+NJ$$ ~ Jump the pointer +n characters from the current position. The 
movement is not stopped when an end of the line is encountered but con- 
tinues to the next line until the number of the required character n is 
reached. 


*OJ$$ - Jump to the beginning of the current line (=0J$$). 


*nGtext$$ - Get the nth occurrence of the specified text. The pointer is 
Placed just after the text string. If n = 1, the 1 can be omitted. 


2.4.4.5 Listing of Text 


In order to be sure that you are at the right place,it is advisable 
to print out the text. 


*L$$ - List out the text from the current position of the pointer to the 
end of the line. 
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*/L$$ - Prints out all the text from the current locations of the pointer 
to the end of the text. 


If, however, you want to see several lines only, the type/L$$ and the 
printing will begin. When you want to stop the printing, type0 (CNTRL 0) 
and the printing will cease. It is important to note that listing does not 
change the position of the pointer. 


2.4.4.6 Changing Characters 
This is actually text modifying. 


*nctext$$ - Change n characters from pointer with specified text. The 
pointer is placed just after the change. 


Example: Suppose you have a line ABC=X+Y and you want to change the 
(+) to (-). First bring the pointer to just before the (+) and the type 
*1C-$$. This means change 1 character immediately after the pointer to -. 


*0Ctext$$ - Replace characters from the beginning of the current line up 
to the pointer with specified text. 


2.4.4.7 Inserting Text 


*Otext$$ - Insert specified text just after the pointer. The use of this 
statement is straightforward - bring the pointer to the required position 
and insert the text. The pointer will be just after the insertion. The 
text can be any legitimate character from a space up to several lines. 
Of course,when inserting more than one line do not forget to type in the 
<cr?. 
2.4.4.8 Deleting Text 

It is possible to delete characters and also lines. 
*+nD$$ - Delete +n characters from the current location of pointer. 
*0D$$ - Delete from the beginning of the current line to the pointer. 
*/D$$ - Delete from the pointer to the end of the text. 


*nK$$ - Delete n lines beginning at the pointer and ending at the nth 
<cR><CLFD (kill n lines). 


*OK$$=0D 


*/K$$=/D 
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2.4.4.9 Chain of Operations 
This feature has the capability to perform chain operations, which 
makes the actual text editing much easier. This is best explained with 
an example: } 
*3A$23I3$2CX1$145$I A2$B$/L$$ 
This means - Advance 3 lines 
Jump 23 characters 
Change 2 characters to Xl 
Jump 14 characters 
Insert A 2 (the A here means one space) 
Back to the beginning of the program 


List the entire program 


Each statement is separated by $ (ESC/SEL) and only when $$ occurs, is 
the chain executed. 


From experience, it was found that it is very easy to err, so it is 
advisable to print before and after the text modification in order to make 
sure that the new text is the desired one. 


2.4.5 Running a FORTRAN Program 
In order to run a FORTRAN program you must have a FORTRAN file, i.e., 


a file written in the FORTRAN language and having the extension .FOR. 
For example,DEMO.FOR. 
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This is your source progran,. 


First, you have to compile your program using the FORTRAN mode, and i 
then you have to link it using the LINK mode and finally run it. The 
whole procedure looks like: 
| 


-R FORTRAN ¢cR) 

*DEMO ,TT:=DEMO ¢ CR D> 

List of program including 

line no. assigned by the compiler 
error messages if any, etc. 

FAC 

-R LINK {CRD 

*DEMO=DEMO/P €CR > 

*AC 

-R DEMO <CR? 


Explanation 


You start in’ KMON as signified by the (.} then you transfer to FORTRAN 
mode and receive the (*). Then essentially you tell the computer to create 
(by default) a file named DEMO.OBJ which will be the compiled version of 
the file DEMO.FOR (also by default). At the same time, you ask for a 
print-out by the TT of the source file with the line numbers given by the 
compiler and all the other statistics. 


If you want to stop the printing at any time, then type~™0. The print- 
ing action will stop but the compilation will continue. When compilation 
is finished, the (*) will be printed. If no error messages occurred then 
you can proceed. First you return to KMON by“C and you get the (.), then 
you transfer to LINK mode. You then ask the computer to create a file by 
the name DEMO.SAV. This is done by the linking program. In this way, you 
get a set of machine instructions from the compiled version. If no error Soe 
messages occurred during linking, then return to KMON, get the (.), and then 
run the program. 
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Actually, R DEMO means by default run the DEMO.SAV. Now, whenever you 
want to run the program just type: 


-R DEMO <CR> 

If you do not want to have a printout in the compilation stage (probably 
to save time and because you are sure of the changes you have done in the 
text), you can type in the FORTRAN mode. 

*DEMO=DEMO <CR > 

This is the same action as before except for the printing. If an error 
message occurs in the compilation or printing processes - correct it 
accordingly. 

Usually, however, a program consists of a main program and several 
subroutines. So a procedure for running this combination is needed. For 
example: 

We have'a file MAIN.FOR which is the main program and two subroutines 
(also as files) SUB1.FOR and SUB2.FOR. The sequence of instruction is as 
follows (assuming no error messages) 

-R FORTRAN (CR > 

*MAIN=MAIN (CR > 

*SUB1=SUB1 <¢CR 7 

*SUB2=SUB2 {CR > 

x “AC 

-R LINK<CR> 

*MAIN=MAIN, SUB1 , SUB2/F (CR ? 
* “co 

«R MAIN <cry 


Here of course each file was compiled separately without printing and 
then they were linked. 
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It is possible, if wished, to compile them together by the statement: 
*MAIN=MAIN, SUB1,SUB2 <CR> 
This creates a file - MAIN.OBJ out of the compilation of the three other 
files, and only it has to be linked. 

It should be noted that it was not necessary to name the .OBJ file 
by the name MAIN but only optional for the sake of convenience. The 
same rule applies in the LINK procedure. 

If printing the program is desired then type 


*MAIN, TT:=MAIN,SUB1,SUB2 <CR> 


This creates a MAIN.OBJ file and gives a printout with the line numbers 
memory assignment, etc. If you want to run the program again, all you 
have to do is type 


-R MAIN <CR? 


If you want to stop execution of a program in the middle of operation 
type twice “™C. 


20 


| 


3. Running of Ray Tracing Program on the PDP 11/40 


As mentioned in the introduction, it is assumed that the reader is 
already familiar with the ray tracing program as described in (1). 
In order to prepare the PDP 11/40 for the ray tracing program, you have 
to type in two statements. 
You start from KMON- 


«SET TTY WIDTH=136 <CR> 
-DAT 08-MAR-78 {CR > 


The first line tells the TT to print lines of 136 characters long | 
(which are needed in our case). The second line is the setting of the 

date. A nine character string starting with the day no., then a minus sign, 
the first 3 letters of the month, a minus sign, and the last two no. of the 
year. In this example, March 10, 1978. 

The date is needed for the program. 

If you want, you can set the internal clock by typing: 

-TIM 09:35 <CR? 


In this example, you set the clock to 9:35 AM. 
For 3:24 PM type 15:24. 


3.1 General Description 


The complete ray tracing program with all the options contains some 
fifty subroutines. Not all of them were transferred to the PDP 11/40. 


List of the subroutines that were chosen together with the reasons 
is given below: 


a. Main body of program includes subreutines 
NITIAL, READW, TRACE, BACKUP, REACH 
POLCAR, PRINTR, RKAM, HAMLTN. 


These were chosen because they are a must, i.e., they are the skeleton 
on which the various options of ionospheric models are added. 


Missing are the subroutines dealing with Plotting since we do not 
have these facilities. 


b. Version of refractive index are 


AHWFWC, BOQWFWC 
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These two were included since the proposed system has to function with real 
ionosphere that includes magnetic field and collisions. The Appleton 
Hartree formula (AHWFWC) is considered as adequate. The Booker Quartic 
version (BQWFWC) is included as a backup although its operation on a sample 
case was not always satisfactory due to machine limitation (as will be 
explained later). j 


c. Versions of electron density models are 


a ener nnnteneree none na 


CHAPX, TABLEX+GAUSEL, TABLE. 


CHAPX was included as a part of the sample case (the one in the 
original version (1)). The sample case as a whole was included in order to 
facilitate a quick check on the whole program in case of doubt. 


Having a Digisonde enables us to measure in realtime the ionospheric 
conditions. The output of the Digisonde is a table of electron density 
values vs. height. So, we need subroutine TABLEX (that goes with GAUSEL). 


Subroutine TABLE was written especially for our purpose. It performs 
essentially the same function as TABLEX and GAUSEL, but it is more 
adapted to our system. 


The interpolation law between the points of height in the table of 
electron density values is parabolic-logarithmic - the same law that is 
used in the Digisonde's signal processing. 


£0 for normal use TABLE should be chosen (does not need GAUSEL) and 
TABLEX+GAUSEL Serve as a backup. 


da. Versions of electron density perturbation models are ELECT1, 
WAVE, TROUGH. 


ELECT1 will serve for the non-perturbation case if so desired. 


WAVE can describe a gravity wave irregularity traveling from north 
to south and is also used in the sample case. 


TROUGH was chosen because it seems to be the most suitable version 
to include tilt of ionospheric layers. 


At the time of writing of this report, it is not known whether this 
subroutine will be adequate or not. The reason is that the transformation 
between the drift attachment measurement (in the Digisonde) and the 


: ionospheric layers tilts has not yet been completed. So, it is possible . 
that writing a new short subroutine to serve as a perturbation model will 
: be needed. 
; 
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e. Version of magnetic field are: 
DIPOLY, HARMNY ’ 
DIPOLY was chosen as a part of the sample case. 


HARMNY (in the original version HARMONY) was chosen as the most complete 
representation of the earth's magnetic field. 


f. Versions of collision frequency chosen are: 
EXPZ and EXPZ2 
EXPZ2 is a part of the sample case. 


Listing of all these subroutines appear in APPENDICES 1-22. Each 
subroutine is written on a separate file having the same name with 
extension, i.e., subroutine READW is in a file named READW.FOR and so 
forth. 


3.2 Main Changes from the Original Version 


In the adaptation process of this program from the original version 
to the PDP 11/40 machine many changes had to be done. The major ones 
are mentioned here: 


- Variable names are allowed to have up to six characters - so in 
case of a longer name in the original value the name was shortened. For 
example, if the original name was HARMONY it was changed to HARMNY, etc. 


- Assignment of Alpha-Numeric string of characters as a variabie's 
value (for printing purposes) is possible on this machine only by the use 
of DATA statement - in contrast to the original version. Furthermore, 
that string of characters should not be greater than four. Therefore, 
longer string had to be divided - like a string of eight characters was 
divided into two strings of four characters each. 


- This machine does not recognize the ENTRY point statement which 
appears in the original version. 


- Use of the BLOCK DATA statement was required in order to assign 
initial values to variables that appear in a COMMON BLOCK, 


- The FORTRAN library of this machine does not contain the ACOS (xX) 
function. Instead, the ATAN(X) was used where needed. 


- Since our hardware does not include either plotter or card punch 
machines, these options were omitted. 
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- The original printout contained azimuth deviation of the ray from 
the original azimuth of transmission. 


For the proposed system real azimuth is more suitable - so it was 
printed instead. 


All the above mentioned changes and more can be seen by comparing 


the listing of the subroutines that are on the DK (see APPENDICES 1-22, 
except #15) to the original version in (1). 


3.3 Running Sequence 


The running sequence is composed of compiling the proper subroutines, 
linking them and then running. The input data should already be in a 
special file. This file will be explained later. 


As an example, we will take the sample case itself and see how to 
run it on the PDP 11/40. We need therefore the following files: 


NITIAL.FOR, READW.FOR, TRACE.FOR, BACKUP.FOR,REACH.FOR, POLCAR.FOR, 
PRINTR.FOR, RKAM.FOR, HAMLTN.FOR. 


These files comprise the basic program. To them, the following models 
have to be added - AHWFWC.FOR, CHAPX.FOR, WAVE.FOR, DIPOLY.FOR and EXPZ2. 
FOR. 
The running sequence is as follows: 
-R FORTRAN {CR > 
*PART1=NITIAL, READW, TRACE, BACKUP, REACH (CR > 
*PART2=POLCAR, PRINTR, RKAM, HAMLTN (CR > 
*PART3=AHWFWC,CHAPX,WAVE CR > 


*PART 4=DIPOLY,EXPZ2 <cr> 


*A 


-R LINK¢ CRD 

*RAYTRC=PART1, PART2 , PART3 , PART4/F CCR > 
#7 

-R RAYTRC (CR > 


This sequence is an extension of similar ones shown in Chapter 2. 
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The division of the subroutines into the different parts is not obligatory, 
and you can change it as you wish. However, you cannot compile together 
all the subroutines (the computer prints an error message), and you cannot 
link together too many subroutines (depending on their size of course). 
This division was found to function properly, so it was adopted. 


After running the program (R RAYTRC), the computer starts execution 
and the following message is typed: 


ENTER DATA FILE NAME* then 


you print the name of the data file and a <cR7 .- The whole message looks 
like - 


ENTER DATA FILE NAME* TEST1.FOR 
By this way, the computer knows where to read data from. The content of 


TEST1L.FOR is shown in APPENDIX 23 (same as the sample case in (1)). 
Printout of this sample case is shown in APPENDIX 27. 


The printout format is explained in Figure 1. 


In APPENDIX 28, printout of a different sample case is shown. Use 
was made of the TABLE subroutine. For this, it was necessary to feed the 
computer with a table of electron density values vs. height. These values 
were computed from the same Chapman model that was used in CHAPX for the | 
sample case. | 
The input data file in this case is TEST3.FOR and is shown in APPENDIX 25. 

The running sequence in this case is as follows: 


eR FORTRAN 


*PART1=NITIAL, READW, TRACE , BACKUP , REACH (CR 7 


*PART2=POLCAR, PRINTR, RKAM, HMLTN (CR 7 


*PART4=DIPOLY , EXPZ2 <cr7 


AC 


«R LINK <CR 7 


q 


Bh *PART3=AHWFWC, TABLE, WAVE <cR 7 
| 
| 
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*RAYTRC=PART1 , PART2 , PARTS , PART4/F <CR » 
AC 

.R RAYTRC <cr) 

ENTER DATA FILE NAME* TEST3.FOR<CR> 


As can be seen, the only differences are in the compilation of PART3 
and in the name of the data file. The results are shown in APPENDIX 28. 
It should be noted that in both cases the running sequence created a file 
PART3.0OBJ. However, there can only be one file by that name on the disk 
The thing that happened is that when the second time occurred, the new 
file was overwritten on the old file by that name. The same rule applied 
to file RAYTRC.SAV created after the linking procedure. If you want to 
keep the old version just call the new version by a different name. 


Two more cases are shown in APPENDICES 29 and 30. In APPENDIX 29, 
results of the combination of the following subroutines are shown 


AHWFWC, CHAPX, WAVE, HARMNY, EXP2Z2. 


i.e., instead of a magnetic dipol representation of the earth magnetic 
field (DIPOLY) a more complete model was used. 


The input data file is TEST2.FOR, shown in APPENDIX 24. 
In APPENDIX 30, results are shown of the following combination: 
AHWFWC, TABLE, WAVE, HARMNY, EXPZ2 


This case differs from the last one by the fact that the electron density 
Model is given by a form of table. 


The input data file is shown in APPENDIX 26. This is the most 
complicated case for it contains TABLE and HARMNY; and it was important to 
verify that it works properly. 


Another case tested was similar to the one shown in APPENDIX 28, but 


instead of subroutine TABLE, subroutines TABLEX and GAUSEL were used. The 
results were similar and therefore, not shown in the report. 
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3.4 Input Data 


As can be seen from APPENDICES 23, 24, 25, and 26, the input data 
file consists of from 1 - 3 main parts. In every case, the first part 
is the one that appears in the original program (1). For convenience 
the rules are copied here: 


The data is entered into a linear array W in the program. 
Each piece of data is stored on one card (one line in the file). The 
first three columns of the card identify the data. Table 2 defines the 
identifying numbers that are subscripts for the linear array - W. The 
. last 56 columns of the card are available for comments. 


The computer program needs angles in radians, whereas people usually 
use angles in degrees. The program is set up for angles in radians but 
putting a "1" in column 18 allows the user to enter the angle in degrees, 
and have the program make the conversion. 


A "1" in column 19 allows the user to enter central earth angles as j 
the great circle distance along the ground in km (the program will calcu- 
late the latitude of a transmitter which is 500 km north of the equator, 
for instance), 


the program expects distances in km, A "1" in column 20 indicates a 
distance in feet. 


For clarity, the read statement and format will be given here (from 
subroutine READW) 


READ (10, 1100) NW, W (NW), DEG, KM, NM, FEET 

1100 FORMAT(I3, El14.7, 411) 
NW - is the identifying number or the subscript. 
W(NW) - is the data itself 
DEG, KM, NM, FEET - are the conversion factors. 
As mentioned earlier, Table 2 explains what to put in. Some values are 
already initiated in the program and do not have to be entered. However, 
if different values are required or desired, then the new values entered 
and read will override the original values that are in the program. 


W 72 and W 81 - 88 do not have to be entered since we do not have punching 
and plotting options in the version. 


ee ee mprieee 
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TABLE 2 


Description of the Input Data for the W Array 


wl = 1. for ordinary ray 
=-l1. for extraordinary ray | | 
w2* Radius of the earth in km 
W3 Height of transmitter above the earth in km | 
w4 North geographic latitude of the transmitter . 
ws East geographic longitude of the transmitter 
W7 Initial frequency in MHz 
ws Final frequency in MHz 
wo Step in frequency in MHz (zero for a fixed frequency) 
Wil Initial azimuth angle of transmission 
wW12 Final azimuth angle of transmission 
wW13 Step in azimuth angle of transmission (zero for a fixed azimuth) 
w15 Initial elevation angle of transmission 
w16 Final elevation angle of transmission 
w17 Step in elevation angle of transmission (zero for a fixed 
elevation) 
W20 Receiver height above the earth in km 
w21 Nonzero to skip to the next frequency after the ray has | 
penetrated the ionosphere ie 
W22* Maximum number of hops | 
W23* Maximum number of steps 
w24* North geographic latitude of the north geomagnetic pole 
w25* East geographic longitude of north geomagnetic pole 
w41* =l1. for Runge-Kutta integration 


=2. for Adams-Moulton integration without error checking 
=3. for Adams-Moulton integration with relative error check 
=4. for Adams-Moulton integration with absolute error check 


W42* Maximum allowable single step error 

w43* Ratio of maximum single step error to minimum single step error 
w44* Initial integration step size in km (step in group path) 
w45* Maximum step length in km 

w46* Minimum step length in km 

w47* Factor by which to increase or decrease step length 

w57 =1. to integrate, =2. to integrate and print phase path 
w58 =1. to integrate, =2. to integrate and print absorption 
w59 =1. to integrate, =2. to integrate and print doppler shift 
w60 =1. to integrate, =2. to integrate and print path length 
w71 Number of steps between periodic printout 

w8l =0. to not plot ray path 


=1. to plot projection of ray path on a vertical plane 
: =2. to plot projection of ray path on the ground 
i W82-88 Parameters used when plotting 
w100-149 Parameters for analytic electron density models 
W150-199 Parameters for perturbations to electron density models 
W200-249 Parameters for analytic magnetic field models 
W250-299 Parameters for analytic collision frequency models 


*These values have been initialized in the main program but may be reset 
by reading them in. See Appendix 1 for initial values. 
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For further details, see APPENDIX 23, note that the last card (line) 
is a blank one, to signify the end of this part. 


APPENDIX 24 shows the input data file for the case of harmonic 
representation of the earth's magnetic field (subroutine HARMNY). 


The portion after the blank card is the input data for HARMNY. This 
data is fixed so it should be left untouched (more details in (1)). 


APPENDIX 25 shows the input data file in case there is a table of 
electron density values and a dipole magnetic field representation (the 
dipole parameters are read in the first part). For more details see 
below. 


APPENDIX 26 shows an input data file similar to the one that will 
probably be used, for it contains the input for the HARMNY and TABLE 
subroutines. 


The first part is the one described in TABLE 2. The second part is 
the fixed data for the magnetic field. The third part contains data of 


electron density values. The second and third parts should not be inter- 
changed. The first line in the third part contains the number of points 


to be read (31 in Appendix 26), the other line contains the height in km 
and electron density (in electrons/cm3). The format is given here (taken 
from TABLE) : 

READ (10, 1000) NOP, ( (HTAB(I), NTAB (I)), I = 1, NOP) 
1000 FORMAT (14/(F8.3.E12.4) ) 
Where: 

NOP - Number of Points 

HTAB - Height of Table in km 

‘NTAB - Electron density (N) of Table in elec/cm>. 
In summary, in order to prepare an input data file, you may copy of of 


the files in APPENDICES 23-26 (according to your needs) and then modify 
it. This way a lot of typing can be saved. 
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3.5 Explanation of Output Format 


When the ray tracing program starts to run, the first printout is 
a reproduction of the input data as read by the computer. This is done 
Samples of this can be seen 


to facilitate easy checks on the input data. 
in APPENDICES 27, 28, 29, and 30. 


APPENDIX 27 is the output of the sample 
TEST1.FOR 

APPENDIX 28 is the output of the sample 
TEST3.FOR 

APPENDIX 29 is the output of the sample 
TEST2.FOR 

APPENDIX 30 is the output of the sample 
TEST5.FOR 


case 


case 


case 


case 


the input to it was 
the input to it was 
the input to it was 


the input to it was 


In the next step, the real calculation begins and a printout of the ray 
trajectory appears. A sample of this is shown in Figure 1 which is self- 
explanatory. This differs from the similar figure in (1) (p. 171) in 


one aspect. 


Instead of printing the: "Azimuth angle of the direction of trans- 
mission in degrees clockwise from the great circle between transmitter 
and ray point" as done in the original version,"Azimuth of ray point from 
transmitter" was printed and the header was changed accordingly. The 
reason for the change was more compatibility with the proposed system. 


The following is a short explanation of the computer messages on 


special points on the ray's trajectory. 


XMTR - at the transmitter 


MIN DIS - at the closest approach point (minimum distance). 


A little elaboration is required to explain the meaning of this - closest 
approach point is defined (in the program) as the point where the wave 
normal is horizontal. It approximates an apogee if it (the min-distance 
point) is above the receiver's height and approximates a perigee if it is 
below the receiver's height. This is only an approximation because this 
condition is found using the wave normal direction which does not necessa- 
rily coincide with the ray's direction in the presence of magnetic field. 


The situation is clarified in Figure 2. 


Whenever the ray's point is one of the two described in the two left 
cases in the figure (and the wave normal is horizontal), the program prints 
MIN DIS. At the two right cases in Figure 2, the program prints: 


APOGEE when at apogee point 


PERIGEE when at perigee point 
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Figure 2: Explanation of the Minimum 
Distance Concept 


RCVR at the receiver's height 
GRND REF ground reflection 
STP MAX exceeded the maximum number of integration 


steps allowed 


PENETRAT ray penetrated the ionosphere 
ENTR ION ray entered the ionosphere 
EXIT ION ray exited out of the ionosphere 
' 
F 
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WAVE REV means wave reversal. This also needs further elaboration. 
The wave normal K is represented in the program in a spherical coordinate 
system (r, ®, 6). At each point on the ray's path, the wave normal 
components K,, Kg, Kg, are calculated. ~ 


The new value is compared to the old one and whenever one of the 
following three conditions occur: 


(Ky) oa ° (Ky) new <o 
(Kg) o14 i (Kg) new <0 
oa * “pnew <° 
the wave reversal message is printed out. 
MAX LAT stands for maximum latitude. 
MAX LONG’ stands for maximum longitude. 
The variables being integrated in this program are r, 8, Gs Ky, Ke: 
Kg, and the independent variable is the group path T. Explicitly two of 


them namely a°/aT ana a5/at ,are integrated. 


At each step of the numerical integration the new value is compared 
to the old one and when: 


C) c?) & ; ; 
(da /AT) sia - (da /8T) | ow <0 MAX LAT is printed 


& D e : : 
(d /8T) 01 a: (a*/aT) aad <0 MAX LONG is printed 


4. Final Remarks 


Comparing the results of calculation as they appear in APPENDIX 27 to 
the original, values show that generally they are very similar. However, 
there are some differences; these occur because of the differences 
between the CDC-3800, for which the program was intially written, and the 
PDP-11/40. 


The main difference stems from the fact that a real variable is 
represented on the PDP-11/40 by 32 bits; whereas it is represented in the 
cpc by 64 bits or so. This means that on the PDP-11/40 there can be 6 or 
7 significant decimal digits compared to 11 or 12 on the CDC, So rounding 
errors are larger in this case. Notice, however, that on the special 
points on the ray's trajectory the differences are small. 


The results between the special points may differ significantly 
because the intermediate printing is done every fixed number of steps 
(S in our case). Step size may differ between the two computers which 
actually print from two different points. 


The special points however, are very similar and the differences stem 
from the rounding errors on the PDP-11/40. 


Two special points stand out as very different-the point where the ray 
enters the ionosphere and the point where the ray leaves the ionosphere. 
The difference stems from the way the logical variable SPACE is calculated 
in subroutines AHWFWC in the two programs. This logical variable is TRUE 
when the ray is in free space and accordingly, the calculation is carried 
out. 


In the original version we have: 

SPACE=REAL (N2) .EQ.1..and.ABS (AIMG(N2)).LT.ABSLIM 
where N2 is the square of the index of refraction and ABSLIM=1.E-5. So in 
simple words, the program decides that the ray is in free space if, at the 
same time, the real part of N2 is equal to one and the imaginary pare of 
N2 is smaller than 1.E-5. 

On the PDP-11/40 this way did not work so it was written aifferently. 

SPACE=ABS (REAL (N2) -1.) LT.0. 5E~7.and. ABS (AIMAG (N2) .LT.ABSLIM 


if instead 0.5E-7.0 was written, then the two statements would have been, 
mathematically speaking, identical. 
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Again,this was necessary because of rounding errors. (Recall that in 
any calculations which include addition and subtraction, the smallest 
meaningful number on the PDP-11/40 is somewhere around E-7.) Even 
though the calculation is done differently in free space (connecting 
straight lines) than in the ionosphere, the actual results, although 
differing at the beginning, quickly converge. This is so because at these 
heights the deviation of the ray from the straight line is negligible. 


Some more statements were changed in the same spirit especially in 
subroutine PRINTR. 


Another example of rounding errors can be seen on the first line of 
the printout of the ray's trajectory. When the ray is at the transmitter 
(XMTR) the range should be zero. However, on the PDP-11/40 it is 
0.0007KM. 


Another kind of limitation encountered, was the reaching of the 
largest number that can be used in the PDP-11/40. This is basically a 
hardware limitation that has to deal with the number of binary bits used 
to represent a real number. 


This problem was first encountered while using the AHWFWC subroutine 
The program stopped and an error message indicating overflow was printed. 
The problem was solved by changing the representation of the algebraic 
expressions where the problem occurred. 


For example: 
Instead of writing A/D**2, A/D/D was written 
or instead of (B**2-4.*A*C), (B-2*SQRT(A) *SORT (C) ) * (B+2*SORT (A) *SORT(C) ) 


was written, etc. For the sample case, this problem was solved as seen 
in APPENDICES 27-30, and it is hoped that it will not occur again in 
the future. 


The same problem occurred while using the BOWFWC subroutine. In this 
case, the problem appeared first in one expression. This expression 
was changed and then the overflow occurred in another expression (the 
first expression for SCALE in APPENDIX 1l). Many ways were tried to 
write this expression without success. In spite of this the subroutine 
was left in the DK for two reasons. One is that further work on it may 
SOlve the problem and the other is that in real situations the subroutine 
May work. 
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